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Abstract 

We study the questions of the existence and mass of the proposed d'(LJ p = 
00 _ ) dibaryon in the quark-delocalization, color-screening model (QDCSM). The 
transformation between physical and symmetry bases has been extended to the 
cases beyond the SU(2) orbital symmetry. Using parameters fixed by baryon 
properties and NN scattering, we find a mild attraction in the IJ P = 00~ 
channel, but it is not strong enough to form a deeply bound state as proposed 
for the d' state. Nor does the (isospin) 1=2 NA configuration have a deeply 
bound state. These results show that if a narrow dibaryon d' state does exist, it 
must have a more complicated structure. 
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1 Introduction 



Quantum chromodynamics (QCD) has been accepted as the fundamental theory of 
the strong interaction. Understanding the low energy behavior of QCD and the na- 
ture of the strong interactions of matter, however, remains a challenge. Lattice QCD 
has provided numerical results describing quark confinement between two static col- 
orful quarks, a preliminary picture of the QCD vacuum and the internal structure of 
hadrons in addition to a phase transition of strongly interacting matter. Phenomeno- 
logical quark model analyses of hadron spectroscopy have also provided useful physical 
information. However, the color structures available in qq and q 3 systems are limited. 
Multiquark systems involve more complicated color structures which can not be stud- 
ied directly in meson and baryon systems. A simple example is given by three gluon 
exchange JTJ and the three body instanton interaction both of which do not contribute 
within a colorless meson or baryon but do contribute to a multiquark system. There- 
fore multiquark systems are indispensable for the full study of the low energy behavior 
of QCD and the structure of strongly interacting matter. 

This report is limited to dibaryon or q e systems. Since Jaffe predicted the H particle 0, 
the study of dibaryons has waxed and waned. At the end of the 1970 's and beginning 
of the 1980's, many dibaryon states were predicted based on the MIT bag model 
and some of them were even claimed to have been observed experimentally. However 
further measurement has almost completely dismissed all of them. Most quark models 
naturally find dibaryon states, except the recent model proposed by Glozman et al [£| , 
which has none|J. 

Experimental signals have been scarce, at best. The search for the H particle has 
been continued for more than twenty years with no indication of its existence. Lomon 
predicted a high mass NN 1=1 1 5'o resonance at around 2.7 GeV|| which does seem 
to be supported by SATURNE pp scattering data|7[. The Moscow-Tubingen- Warsaw- 
Uppsala collaboration || has claimed a narrow dibaryon resonance, d', centered at 2.06 
GeV with a small width of 0.5 MeV; the preferred quantum numbers are IJ P = 00~. 
In contrast to all the other cases, the d! dibaryon, with a mass as small as 2.06 GeV 
and IJ P = 00~, is hard to accommodate by the available quark models||. Newer 
experiments using simple systems have not confirmed the existence of a d' signal flTUf . 

A new quark model, the quark-delocalization, color-screening model (QDCSM), has 
been developed with the aim of understanding the well known similarities between 
nuclear and molecular forces despite the obvious energy and length scale differences [ pi]] . 
The model predicts a small mass narrow dibaryon resonance d* with IJ P = 03 + , M~2.1 



GeV, T(NN) ~ 1 MeV|L3]]. Although the model has not been applied to the study 
of baryon resonances above the ground state flavor octet and decuplet, where other 
models have shown good accuracy |12|], we have found that it tends to underestimate 
absolute orbital excitation energies. However, since our calculations are structured to 



1 



self-consistently compare dibaryon states to our calculated two baryon thresholds, this 
should not lead to an underestimate of the d' dibaryon mass, within the model. 

This paper reports the results of our study of the d! dibaryon using the QDCSM. 
Sect. 2 gives a brief description of the Hamiltonian and wave function in QDCSM; Sect. 3 
describes the calculation method. For this, a very useful extension of the transformation 
between the physical bases and symmetry bases beyond the SU(2) orbital symmetry is 
presented. The results are presented in Sect. 4. The final section provides a summary. 



2 Model Hamiltonian and wave function 



The details of the QDCSM can be found in Ref.|TT], [Tj|. Here we present only 
model Hamiltonian and wave functions used in the calculation. 

The Hamiltonian for the 3-quark system is the same as the usual potential model, 
the six-quark system, it is assumed to be 



the 
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if i, j occur in the same baryon orbit, 
if i,j occur in different baryon orbits, 



where all the symbols have their usual meaning except the confinement potential 
which will be explained below. 

The wave function of the six-quark system is written as 
where is the 3-quark system wave function, 



and 



^ B2 = [[^(4)^(5)^(6)]S/ 2 s 2 P* ,Xc(456). 
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(2) 

(3) 
(4) 



The single particle orbital wave functions are delocalized, as 



+ e s <j) R )/N s , 
*>r + e s (p L )/N s , 

1 X 3/4 
Tib 2 ) 



(5) 



N, = [l + E ; + 2 Eie *""Y /2 

if the particle is in s-wave and as 

= (4>l + e p (j)R)/N p , 
tpR = {<t>R + e p (f> L )/N p , 



^) 3/4 ^\r + s /2\Y lm e^ +g ^ 2 , (6) 

^) 3/4 ^|r-V2|r lm e-^^ 2 ) 2 , 
rov b 

N s = [l + el + 2e p (l-^ 2 )e^ b2 ] 1/2 

if the particle is in p-wave, where sq is the separation of the two clusters. 

Although we use a potential model language in the description of the confinement 
potential in Eq.(l), our calculations employ an extended effective matrix element 
method. Even though understanding of quark confinement is limited, the use of a two 
body interaction to describe the quark confinement may well be highly oversimplified, 
especially for multiquark systems. For example, the three-gluon and three-body in- 
stanton interactions mentioned above cannot be expressed in terms of two-body 
interactions and the full nonperturbative QCD interaction contains additional varieties, 
such as condensates, which can not be expressed in terms of two body interactions ei- 
ther. Our confinement "potential" as defined in Eq.(l) has the meaning of the usual 
interaction potential only in the asymptotic region where the overlap of the <f>L and 
<j)R wave function orbitals is negligible. In the interaction region, it defines a recipe for 
determining the effective matrix elements of the Hamiltonian of a six quark system. 
(More details will be provided in the next section.) Our goal is to include the major 
portion of nonperturbative many-body quark interactions by this method. 

We refer to the six-quark wave function defined in Eq.(|2|) as the physical basis because it 
has an obvious physical meaning in terms of a pair of baryons. In our channel coupling 
calculation, all of the colorless p-wave excitation baryon-baryon channels compatible 
with the d' quantum numbers IJ P = 00~ are included. Although the six-quark wave 
function, Eq.(^), appears to have the form of q 3 — q 3 clusters, in fact the other clusters 
such as g 6 , q 5 q and q A q 2 have been included because of the quark derealization. Hidden 
color channels have been omitted because we do not know anything about the mass and 
interaction of colorful baryons and because a hidden color channel can be expressed in 
terms of a sum of colorless channels \TA\. 
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3 Calculation method 



The physical basis used in Eq. (|2|) is not convenient for the calculation of matrix elements 
of the six-quark Hamiltonian. In order to simplify the calculation, the physical basis 
(cluster basis) is expanded in terms of symmetry bases first. Next, we use the powerful 
fractional parentage (fp) expansion method to calculate the matrix elements between 
symmetry basis components of six-quark system. Finally, the matrix elements between 
physical basis components are obtained by reversing the transformation between the 
physical and symmetry bases. 

The transformation between the physical basis and symmetry basis has been studied by 
Harvey P^, Chen|T^], and ourselves |T7|] in the case where the orbital symmetry of the 



three-quark cluster was limited to [3] and SU X (2) orbital symmetry was assumed. In the 
case of interest here, there is a p-wave quark with respect to the right- or left-center, in 
addition to right- and left-centered s-wave quarks. The orbital symmetry group needs 
to be extended to SU X (4). (Only one state among the three p-wave states is included 
in the calculation of the transformation coefficients, since the spin-orbit coupling is 
a trivial one.) The orbital symmetry of individual baryons should include both [3] 
and [21] configurations. This requires and extension of the transformation method 
studied previously. Here we develop two new methods to calculate the transformation 
coefficients. 



In the symmetry basis, the group-chain classified basis is denoted by 

W]w x \ 



\®s) 



[a}W c [ii]IM I SLJM_ 



(7) 



where [u], [a] and \p] represent the symmetries of orbital, color and spin-isospin degrees 
of freedom, and W X ,W C are Weyl tableaux for [u] and [a]. Other symbols have their 
usual meanings. The group chain we use here is 

SU XCTa {36) D SU X {3) x SU CTa (l2) D SU C (3) x SU Ta (A) D SU T {2) x SU a {2). (8) 

The physical basis is constructed from two 3q states, 

\*p)=A[B 1 B 2 } [ S j Mj} (9) 

where B\ and B 2 represent the wave functions of two three-quark clusters. The physical 
basis defined here is a little different from the ty§ q defined in Eq.(^). The two bases 
are related by Racah coefficients. Fortunately, for the quantum numbers of d', all of 
the coefficients equal unity. From now on, we replace by \l/p. 

The physical basis can be expanded in terms of the symmetry basis as 

l*p) = E c P -s\®s). (io) 

uWxjJ, 
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By using the unitary condition of the symmetry basis, the expansion coefficients can 
be expressed as, 



a 



p—s 



($s|^ P ). 



On the other hand, the six-quark symmetry basis can be expanded into two three-quark 
clusters by application of the fractional parentage technique: 



[i 6 ],MP] 
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where 1, 2 stand for v\, z/ 2 , W^, W X2 , fii, /12, etc. The last line is just the wave func- 
tions of two three-quark clusters. Combining with the relevant Clebsch-Gordan (CG) 
coefficients, we have 



3 ]KPi 
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(13) 



We next apply the inter-cluster antisymmetrization operator 



A=\ — 
20 



pdnter-cluster permutations 



to Eq.fllB]). Because of the total antisymmetrization of the six-quark symmetry basis, 
when the operator A is applied to the left-hand side of Eq. (pT3|) , it produces 20 copies. 
In this way, we obtain 



[u]W x 
[a]W c [^} IM i SLJM J 



l°\Mi>\ 



p],WM 



C 



[u]W x 



1,2 



Substituting Eq.(14) into Eq . (|TT|) , and making use of the orthonormal property of the 
physical bases, we obtain the expression for the transformation coefficients, 



(15) 



This expression differs from the transformation coefficients obtained before [16, [L7j. It 



has two "new" factors: C, 



l%[u][u] 

[13] M [/>!], [13] \u 2 ]\u 2 ] 



and Ci 



v\W x 



In the limit of SU X {2) 
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with [ui] = [u 2 ] = [3], the product of these two factors is a constant with the value, 
which returns us to the expression of Ref. |U| [L7] . 

Recently Chen[jlS|| also proposed a method to generalize the transformation between 
physical and symmetry bases. Our new method is based on his expressions in Eqs.(9- 
20a) of Ref.pl, 



A 



[vi]a h 



\v 2 \b h 
(T2H2I2S2 



UfJ, 



v{\ [0-1] [m] , [u 2 ] [ct 2 ] [^2] [m]hSi, [jtta] 



c 



Ms 

[v]a h b h 
afilS 

(16) 

To extend this to the case fo interest here, first set a? 1 to W Xl , b^ 2 to W X2 . The 
symmetry basis used in Eq. (|T^) is not the one which we defined previously. Rather, it is 
an SU4 D SU2 x SU2 irreducible basis, which is a non-standard SU4 basis, with respect 
to the orbital symmetry and must be expanded in terms of the standard Gel'fand bases 
of SU/l and the Yamanouchi bases of 5*6. It is only in the special SU X {2) case that the 
Weyl tableau, W x , is automatically fixed by the orbital symmetry [v\. In general a 
further transformation from non-standard to standard SU„ bases is needed. 



M Nh] \ v 

m W Xl W x2 J 



m,W x 




w xi w x2 



(17) 



where the transformation coefficients, termed subduction coefficients (SDC) of ST/4, are 
independent of m. From the SDC and isoscalar factor tables of SU4 and CG coefficient 
tables of SU 2 , it is easy to verify that the SDC used here is just the product of v20, 



C 



and C, 



'v\W x 



HHpj.iisiw^j ^[ V1 ]W X1 ,W 2 ]W X2 - So the two methods give the same transformation 
coefficients as they should. 

For the quantum numbers of d', there are six physical channels, (taking only color 
singlet baryons into account,) which are given in the first column of Table I, where: 
N denotes the state of nucleon with quantum numbers [i/J = [3], ij = |,Sj 



0,ji = 21 A denotes [u^ = [3],/, 
I q. - I 1. - 1 a. — I. 

2 , u, 2 J l « x ) Ji 2 ' 



3 C. _ 3 7 

2 5 2 2 ' 2 

A^ 2 * denotes [z/J = [21], J, 



1 I = 

2 ' 1 

0, ^ = |; N* denotes [v t ] = [21], = 

21 = 2>h = ^-iji = 2' all< ^ ^* 



1, Ji 

denotes [z/j] = [21], /j = |,Si = = l,jj = |. The corresponding symmetry bases 
are given in the first row of Table I. There are thirteen bases. The reason we have a 
6 x 13 table is that the hidden color channels have been omitted. The d' is assumed 
to be a linear combination of these six physical channels. The combination coefficients 
are determined by diagonalizing the Hamiltonian in the 6-dimensional space. 

To calculate the matrix elements of the six-quark Hamiltonian between symmetry 
bases, we use a fractional parentage expansion. Then only the two-body matrix ele- 
ments and four-body overlaps are required. Details can be found in Ref. ||17||. By using 
the transformation coefficients derived above, the matrix elements of the six-quark 
Hamiltonian between physical bases can be obtained straightforwardly. 
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In the calculation of the two body matrix elements of the confinement interaction, 
we assume taht the normal confinement interaction should be used for (LL\V C \LL) , 
(RR\V \RR), (LR\V C \LR) , while the color-screened confining interaction should be 
used for (LR\V \RL) and all others. Here 

(LL\V C \LL) = (Mfi)Mr2)\V£ 2 \Mri)Mr2)). 

This recipe fixes the matrix elements of the confinement interaction and distinguishes 
our extended effective matrix element approach from a two-body potential-model ap- 
proach. It reduces to the usual two-body potential-model for a single hadron and there- 
fore maintains the successes of the constituent quark model for hadron spectroscopy. It 
also reproduces qualitatively correct phase shifts in ten channels: NN ST = 10, 01, 00, 11, 
iVA ST = 0|, l|, NT, ST = 0~, 1~, 0|, l|, with only one additional adjustable param- 



eter, the color screening constant //[ll, 20]. It is the only model which demonstrates 



a similarity between nuclear and molecular forces. We take these successes as an indi- 
cation that this effective matrix element approach includes significant components of 
correct physics. 

For a preliminary study, we use the adiabatic approximation. That is, for each separa- 
tion So, we determine the energy of the six-quark system by the variational condition 

dE 

—^ = 0, l = s,p (18) 
oei 

The effective potential between two baryons is obtained by a subtraction [[2"T|], 

V e (s )=E 6 (s )-E 6 (oo). (19) 

and the mass of the six-quark system is estimated by the formula, 

M 6 = m 1 +m 2 + V e + E (20) 

where mi,m2 are the masses of the two baryons in isolation and Eq is the zero-point 
energy of the pair, Eq = 3^7, with m the reduced mass of the two baryons. 



4 Results and discussion 



The model parameters, which are fixed by the baryon spectrum and NN scattering, 
are: 

m u = m d = 313 MeV, b = 0.603 fm, a = 25.13 MeV/fm 2 , a s = 1.54, // = 1.6 fm~ 2 . 

Both single channel and channel coupling calculations were carried out. The results 
are shown in Table II. From these results, we observe a mildly attractive interaction 
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between N and N*, ~ 100 MeV. Since the delocalization is approximately 1, this is a 
six-quark state, not a two-baryon state. However this attraction is not enough to form 
a deeply bound state such as d! . Channel coupling adds a little more attraction (~ 10 
MeV), does not change this conclusion. 

These results are expected in the QDCSM. From our previous study of the effective 
potential between baryons [EI], we found that as a general trend there are strong 



attractions in decuplet-decuplet channels, but the attraction in octet-octet channels is 
weak. The excited N* state of the nucleon is still in a flavor octet, so the attraction 
between N and N* would not be expected to be strong, although the existence of the 
p-wave quark adds a little more attraction. The mass of AA* channel is much larger 
than that of the NN* channels, so the effect of channel-coupling to it is small. 

The effect of the harmonic oscillator parameter, b, of the quark orbital wave function 



has also been studied. G. Wagner et al. [22| concluded that it is impossible to describe 



the dibaryon, d', with the same parameters as those used for single baryons. In par- 
ticular, they found it impossible to get a d' mass as low as 2.06 GeV if the same b is 
used for both 3-quark and six-quark systems. To check their statement, the results 
obtained with a larger value of b in the six-quark calculation are also shown in Table 
II. By varying b, the attraction between N and N* could be increased. However, if a 
different b is used in the 3-quark and six-quark calculations, the subtraction procedure, 
shown in Eq . (|i~9"D to obtain the effective potential, is no longer reliable. A dynamic 
calculation is needed. It is possible that a deeply bound d! state might be obtained 
in the QDCSM if a larger parameter b is used for the six quark system. However, the 
physical meaning of such a model calculation is obscure. 

The possibility that the d! has isospin 1 = 2 has been discussed in the literature. 
We have carried out both adiabatic and dynamical IJ P = 20~ NA single channel 
calculations. With the same model parameters, the adiabatic calculation produces a 
mass of 2173 MeV for the d' . A dynamical calculation with eight Gaussians (spanning 
cluster separation coordinate values from 0.6 to 4.8 fm) produces a mass of 2191 MeV 
and one with fifteen Gaussian (spanning cluster separation coordinate values from 0.6 
to 9.0 fm) produces a mass of 2178 MeV. These trends strongly suggest that this is a 
scattering state, not a bound state. So we conclude that no state with a mass as low 
as 2.06 GeV can be obtained in the NA channel in the QDCSM. 



5 Summary 

We have studied the proposed d' state in the QDCSM using the physical bases NN* 
and AA*. Our results show that a mild attraction does develop in the IJ P = 00~ state, 
but this attraction is not strong enough to form a deeply bound state with a mass as low 
as 2.06 GeV. All of the p-wave excitation baryon-baryon channels with (i'-compatible 
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quantum numbers have been included in our channel coupling calculations. Because 
the QDCSM includes quark derealization, all of the configurations, q 6 , q 5 q, q 4 q 2 and 
q 3 q 3 , are included in our model space. 

Nor do we find such a low mass state in the 1 = 2 NA channel in our model approach. 

Including the d' results of the Faessler group|| |2^], we conclude that, if the d! is 
experimentally verified as a real dibaryon resonance, it must have a more complicated 
structure than studied here, such as including q 7 q components. 

Special thanks are due to Jin-Quan Chen for constructive discussions. This research 
is supported by the National Science Foundation of China, Fok Yingdung Educational 
fund, National Science Foundation of Jiangsu Province and in part by the U.S. De- 
partment of Energy under contract W-7405-ENG-36. 
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Table 1: The transformation coefficients between physical and symmetry bases. 
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Table 2: Effective potentials and masses of six-quark systems. 





fe(fm) 


s (fin) 


e s 




£ 6 (s )(MeV) 


£ 6 (oo)(MeV) 


K(MeV) 


M 6 (MeV) 


s.c. 


0.603 


1.0 


1.0 


1.0 


2376 


2473 


-97 


2466 


c.c. 


0.603 


1.0 


0.6 


1.0 


2364 


2473 


-109 


2454 


s.c. 


0.85 


0.9 


1.0 


1.0 


2123 








c.c. 


0.85 


0.9 


1.0 


1.0 


2121 
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